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Abstract. We report results on quarkonium production from the STAR experiment at the Relativistic 
Heavy-Ion Collider (RHIC). J/ip spectra in p+p and Cu+Cu collisions at ^/s NN = 200 GeV with transverse 
momenta in the range of 0.5-14 GeV/c and 5-8 GeV/c, respectively, are presented. We find that for pr > 5 
GeV/c yields in p+p collisions are consistent with those in minimum-bias Cu+Cu collisions scaled with 
the respective number of binary nucleon-nucleon collisions. In this range the nuclear modification factor, 
Raa, is measured to be 0.9 ±0.2 (stat). For the first time at RHIC, high-pr J/^-hadron correlations were 
studied in p+p collisions. Implications from our measurements on J/ip production mechanisms, constraints 
on open bottom yields, and J/tj) dissociation mechanisms at high-pT are discussed. In addition, we give 
a brief status of measurements of T production in p+p and Au+Au collisions and present projections 
of future quarkonia measurements based on an upgrades to the STAR detector and increased luminosity 
achieved through stochastic cooling of RHIC. 
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The dissociation of quarkonia due to color-screening of 
their constituent quarks in a Quark-Gluon Plasma (QGP) 
is a classic signature for deconfinement in relativistic heavy- 
ion collisions pQ . The suppression of the various J/ip and T 
states is determined by their binding energy and the tem- 
perature in the plasma. Results from the PHENIX experi- 
ment at RHIC show that the suppression of J/tp as a func- 
tion of centrality (the number of participants) is similar 
to that observed by NA50 and NA60 at the CERN-SPS, 
even though the temperature and energy density reached 
in these collisions is significantly lower than at RHIC [H 
[2]. This indicates that at RHIC energies additional mech- 
anisms countering the suppression, such as recombination 
of charm quarks in the later stage of the collision, may 
play an important role; they will need to be studied sys- 
tematically before conclusions from the observed suppres- 
sion pattern can be drawn. In most theoretical models, 
the dissociation of quarkonia states is computed at rest 
relative to the QGP p] . Recently, techniques based on the 
AdS/CFT duality have been utilized to study the dissocia- 
tion of quark-antiquark pairs with large velocities relative 
to the strongly coupled QGP. Calculations in this frame- 
work show that bound states of heavy fermion pairs (an 
analog of quarkonium in QCD) have an effective disso- 
ciation temperature that decreases with px as l/y^y [3], 
implying an increasing suppression, i.e. decreasing Raa, 
with pt ■ To test this conjecture measurements of J/ip 
Raa above pr > 5 GeV/c are needed where the effective 
J/ip dissociation temperature is expected to be lower than 



website: http://www.star.bnl.gov 



the temperatures reached in RHIC collisions (~1.5 T c ) [3J 

urn- 

Understanding J/ip production requires knowing what 
what fraction of J ftps are produced from (i) gluon and 
heavy-quark fragmentation, from (ii) decay feed-down of 
B mesons and \c states, and (in) what fraction origi- 
nates from direct production either through a color-octet 
or color-singlet state. Our understanding of J/ip produc- 
tion mechanisms has gone through several cycles in recent 
decades [S1[SJ[7]. There is no convincing model explaining 
all the major features of the existing data in e + e~ an- 
nihilation and hadron-hadron collisions. The color-singlet 
model (CSM) [5j[6] underpredicts the J/ip spectra by an 
order of magnitude in p + p collisions at y/s = 1.8 TeV [5]. 
The color octet (COM) and color evaporation model (CEM) 
were proposed to explain the production yields, but fail 
to explain the recent measurements of large spin align- 
ment (polarization) from the same experiment [5,2. In 
addition, J/ip is found to be dominantly associated with 
open-charm pair production in e + e~ annihilation at s/s 
= 10.6 GeV [T. These later findings suggest that yields 
calculated from leading-order pQCD may not be the dom- 
inant contribution to the J/ip production. UA1 and DO 
have analyzed the J/i/i-hadron correlations in p + p colli- 
sions to separate J/ip sources at high pt [8] from \c decay 
without a near-side correlation and from B hadron decay 
with a strong near-side correlation [8]. 

In this paper we report on the measurement of J/ip 
production at midrapidity with the STAR experiment in 
p+p and Cu+Cu collisions at ^/s^ = 200 GeV with trans- 
verse momenta in the range of 0.5-14 GeV/c and 5-8 GeV/c, 
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Fig. 1. The e e~ invariant mass distribution in p+p collisions 
and Cu+Cu collisions at ^/s NN = 200 GeV. The solid and 
dashed lines represent the distribution of un-like and like-sign 
pairs, respectively. 



respectively. In addition, J/V'-hadron correlations, orig- 
inally proposed and studied by UA1 [5] are presented. 
The technique used is analog to that deployed in hadron- 
hadron correlations studies in STAR [9] . We will also briefly 
report on the status of measurements of the T in p+p and 
Au+Au collisions and present projections of yields that 
will become available after completion of the STAR de- 
tector upgrades and the increase of RHIC luminosity that 
will be achieved once stochastic cooling for ion beams is 
implemented. 

In the analyses reported here, the J / 1\) and T are recon- 
structed through their decays into electron pairs, QQ — > 
e + e~. The large acceptance of the tracking system in the 
STAR Time Projection Chamber (TPC) [TO] and the elec- 
tron trigger capability from the Barrel Electromagnetic 
Calorimeter (BEMC) [10] with \rj\ < 1 are very well suited 
for such analyses. At STAR, both the TPC and BEMC can 
provide electron identification [10]. At highly, the BEMC 
is a very powerful tool for electron identification and can 
be used for online triggering to enrich the electron sample. 
At moderate px, the TPC provides sufficient energy- loss 
(dE/dx) resolution to identify electrons efficiently. 

Various online trigger schemes were deployed to max- 
imize the recorded yield. The \ow-px J/tp trigger is based 
on two spatially separated BEMC towers with Et > 1-2 
GeV. A higher level trigger makes a final selection based 
on the approximated invariant mass of the pair. However, 
this trigger can only be efficiently used in p+p collisions; 
in A+A collisions this trigger does not provide sufficient 
discrimination power due to the high hit occupancy in the 
BEMC. The high-pr J/ip trigger requires only one tower 



above a given high-energy threshold. For the data pre- 
sented here the thresholds were: Et > 3.5 GeV (p+p run 
5), E T > 3.75 GeV (Cu+Cu run 5), and E T > 5.4 GeV 
(pp run 6), respectively. The highly efficient T trigger is 
based on a single high-tower signal (Et > 4.5 GeV) and 
a subsequent pair selection with invariant mass cut per- 
formed in an online high-level trigger system using the full 
BEMC tower data. 

For the T measurement the sampled luminosities were 
9 pb _1 in p+p collisions and 12 pb^ 1 (p+p equivalent) in 
Au+Au collisions. The \ow~pt J/ip is from a dataset with 
0.4 pb _1 of p+p luminosity. The data used for the high-p^ 
J/ip analysis was recorded during the p+p and Cu+Cu 
runs in 2005 and the p+p run in 2006 using the trigger 
setup described above in coincidence with a minimum bias 
trigger which required a coincidence between the two Zero 
Degree Calorimeters (ZDCs). The integrated luminosity 
was ~ 2.8 (11.3) pb _1 for p+p collisions collected in year 
2005 (2006), and ~ 860 /zb" 1 (3 pb" 1 p+p equivalent) 
for Cu+Cu collisions. In Cu+Cu data, the most central 
0-60% of the total hadronic cross-section was selected by 
using the uncorrected charged particle multiplicity at mid- 
rapidity (\r]\ < 0.5) P2||3|. 

For the high-py J/ip analysis, we required one high-p^ 
electron identified with the combination of BEMC tower 
energy, shower shape (from shower-max detectors embed- 
ded in the BEMC), and dE/dx measured in the TPC. The 
Pt cut for the second electron was substantially lower and 
only dE/dx in the TPC was used for identification. 

Figure [TJ shows the invariant mass spectra for the high- 
Pt dielectron sample in p+p collisions from runs 5 and 6 
(2005 and 2006) and in Cu+Cu collisions (run 5). The 
combinatorial background, derived from like-sign pairs, is 
depicted by dashed lines. The J/ip signal is extracted in 
the mass window 2.9 < M inv < 3.2 GeV/c 2 . Due to the 
high-pT cuts used in this study the signal-to-background 
ratio (S/B) is exceptionally large. In p+p collisions we 
obtain S/B = 22/2 (40/14) for run 5 (6) and 17/23 in 
Cu+Cu collisions. The pt coverage in p+p and Cu+Cu 
collisions taken in run 5 is 5 < pt < 8 GeV/c. For p+p 
collisions in run 6 the larger integrated luminosity, the 
larger BEMC coverage, and optimized trigger thresholds 
allow us to extend the pt reach to 14 GeV/c. Figure [2] 
shows the fully corrected J/tp invariant cross-section B ee x 
Ed 3 a /dp 3 as a function of pt in p+p collisions for all 
three datasets. Our spectrum is in good agreement with 
results from the PHENIX experiment, which is plotted for 
comparison [11) . 

Figure [3] shows the J / ip nuclear modification factors 
Raa as a function of pt in 0-20% and 0-60% Cu+Cu 
collisions calculated from PHENIX 12J and STAR mea- 
surements. The data suggest that Raa is rising towards 
unity for pt > 5 GeV/c, although the large errors cur- 
rently preclude strong conclusions. A combined fit to all 
high pt data points above 5 GeV/c gives Raa = 0.9 ± 
0.2 (stat), consistent with unity and 2er higher than that 
at 1ow-pt (Raa ~ 0.5). This result is in contradiction 
with expectations from AdS/CFT-based [3] and Two- 
Component- Approach [14] models that predict a decrease 
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Fig. 2. J/ip invariant cross-section times the di-electron 
branching ratio as a function of pr in p+p collisions at y/s 
— 200 GeV. The lines correspond to pQCD with Color Evap- 
oration Model (CEM) calculations [2]. 
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Fig. 4. J/?/>-hadron azimuthal correlations after background 
subtraction. The histogram is J/i/;-hadron from B decay with 
full-event simulation from the PYTHIA event generator. 
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Fig. 3. J /ip Raa as a function of pr- The solid line represents 
the fit to all the data points at 5 < pr < 10 GeV/c yielding 
Raa = 0.9 ± 0.2 (stat). The curves are "Two Component Ap- 
proach" model prediction [14] . The boxes in the right show the 
normalization uncertainty for 0-60% (left) and 0-20% (right) 
collisions. 



in Raa with increasing px- A similar trend was also ob- 
served by NA60 Collaboration in In+In collisions at ^/s NN 
= 17.3 GeV [T5|. There, however, Raa reaches unity at 
considerably smaller px than at RHIC, suggesting that 
the effects are most likely of a different physics origin. Our 
results could indicate that other J/ip production mecha- 
nisms that counter the suppression such as recombination 
and formation time effects }16j play an increasingly dom- 
inant role at higher px . The small suppression of the J / ip 
at high-pT stands in contrast to a substantial suppression 
of open charm production at similar px |17j . 

The large S/B ratio of the high-py J /if) data in p+p 
collisions allows us, for the first time at RHIC, to study 



J/'0-hadron correlations, which can potentially provide 
important constraints on the underlying J/ip production 
mechanisms. Figure [4] shows the azimuthal angle correla- 
tions dN/dAcp per J/ip between J ftps with px > 5 GeV/c 
and charged hadrons wither > 0.5 GeV/c. No significant 
yield for near-side correlations (A<f> ~ 0) is observed. This 
is in stark contrast to the case of hadron-hadron correla- 
tions [3]. 

Monte-Carlo simulations (PYTHIA 6.3.19) 18J show 
a strong near-side correlation dominantly due to the feed- 
down of J/tp from B- meson decays [8], B — ► J/ip + X. 
Thus, the comparison of the measured near-side yields 
with PYTHIA simulations allows one to infer the fraction 
of J/ips originating from B meson decays that contribute 
to the observed yield, albeit in a model dependent way. As- 
suming no other contribution to the near-side, we obtain 
an upper limit of 15% of the £?-meson feed-down to the 
inclusive J/i/j cross-section at px > 5 GeV/c. The shape 
on the near-side is not matched well. A detailed compar- 
ison of the spectral shapes is not currently possible due 
to the limited statistics. This will requires data samples 
which we hope to collect in the upcoming RHIC run. 

It is worth noting that the direct comparison of our 
inclusive high-pr J/ip spectra with the feed-down spec- 
tra from pQCD calculations provides less constraints than 
the method described above. Using NLO predictions for 
B meson production |19j and decay parameter from mea- 
surements by the CLEO collaboration [20] we obtain an 
upper limit for the B meson feed-down fraction of ~ 40% 
for px > 5 GeV/c. This is in large parts due to mass and 
scale uncertainties in the NLO calculations [2]. Further 
studies of J/i/'-hadron correlation and J/ip cross-section 
with higher statistics will allow us to constrain the B 
cross-section substantially in the future. 

Figure [5] shows a projection of Raa for an integrated 
luminosity of 5 nb _1 for Au+Au collisions and 200 pb _1 
for p+p collisions with an upgraded STAR detector. The 
key upgrades are: (i) the full barrel time-of-flight detector 
to improve the electron identification in the low-pr region 
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Fig. 5. Projection of Raa at RHIC II luminosity, in which 
STAR will sample 5 nb -1 and 200 pb~ in Au+Au and p+p col- 
lisions, respectively. The result is compared to Raa at 1ow-pt 
obtained by PHENIX in run 4, corresponding to a combination 
of sampled luminosities of 157 /xb _1 Au+Au collisions and 2.6 
pb~ p+p collisions. 



(px < 3 GeV/c) and (ii) the upgrade of the TPC readout 
electronics in conjunction with an upgrade of the data 
acquisition system which will increase our data-taking rate 
tenfold and eliminate dead-time losses for rare triggers. 
This will allow STAR to make full use of the increased 
luminosity of RHIC once the stochastic beam cooling is in 
place. 

In addition, with the improved electron identification 
capabilities and higher data-taking rate the measurement 
of J ftp elliptic flow becomes feasible. Simulations show 
that with the statistics of 10 9 minimum bias Au+Au events 
we will be able to measure v-i with < 1% statistical error. 
This will provide stringent tests of quark coalescence and 
charm flow [2J. 

Many of the difficulties that plague the interpretation 
of the observed charmonium suppression can be avoided 
when studying bottomonia. Due to the small production 
cross-section, recombination effects become negligible and 
absorption by hadronic co-moving matter is unimportant 
|21j . The small cross-section, however, also makes bot- 
tomonium states extremely difficult to measure. While the 
T(1S) is predicted to be not suppressed at RHIC and even 
LHC energies, the T(2S) (T(3S)) state is assumed to disso- 
ciate at temperatures similar to that of the J/ip (ip'). The 
separation of the three states, while possible in STAR, re- 
quires substantial statistics which will only become avail- 
able with the the RHIC luminosity upgrades. Detailed T 
simulation and projections can be found in [2j. 

The STAR experiment already reported on the first 
RHIC measurement of the T(1S+2S+3S) cross-section at 
mid-rapidity in p+p collisions at s/s — 200 GeV. We find 
BR ■ dcr/dy\ y= o = 91 + 28 (stat.) ± 22 (syst.) pb, which 
is consistent with the world data and NLO pQCD calcu- 
lations in the CEM [52]. The first ever measurements of 
T mesons in Au+Au collisions at y/s^ — 200 GeV are 
underway. We observe a stable signal that will allow us 



to get first information on the nuclear modification factor 
of the T. This will be complemented by measurements in 
d+Au collisions taken in 2008. 

In summary, we report measurements of J/ip spec- 
tra from 200 GeV p+p up to transverse momenta of 14 
GeV/c and from minimum bias Cu+Cu collisions at high 
Pt (5 < px < 14 GeV/c) at mid-rapidity. The J/ip nuclear 
modification factor Raa in Cu+Cu at pr > 5 GeV/c is 
0.9 ± 0.2 (stat) and is about 2a above the values at low 
Pt measured by PHENIX [12] . The study of J / ip -hadron 
azimuthal correlations show an absence of near-side cor- 
relations. Using PYTHIA simulations we derive an upper 
limit for the fraction of feed-down of J/ip mesons from 
B meson decays of 15% at pr > 5 GeV/c. Near future 
upgrades to the STAR detector and the RHIC machine 
will allow us to make substantial contributions to the un- 
derstanding of quarkonia production and provide detailed 
information of their interaction with the hot dense matter 
created at RHIC. 
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